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Magnetic and electric properties of the Hubbard model with binary alloy disorder are studied 
within the dynamical mean-field theory. A paramagnet-ferromagnet phase transition and a Mott- 
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induced enhancement of the Curie temperature is demonstrated and explained by the effects of band 
splitting and subband filling. Different quantum phase transitions driven by changes of the alloy 
concentration are identified. 
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I. INTRODUCTION 

There is a great concern about the nature of itinerant 
ferromagnetism in correlated electron systems with disor- 
der, such as, e.g., doped manganites (Lai-xSr^MnOa)^ 
alloyed ruthenates (SrRui-^Mn^Oa)^ or alloyed Kondo 
insulators (FeSii-^Gea)!^ Moreover, it is of funda- 
mental importance for industrial applications to pre- 
cisely control the Curie temperature of different ferro- 
magnetic alloys or diluted magnetic semiconductors, e.g. 
Gai_ x Mn x As&£ 

In typical transition metals (e.g. Fe, Ni or Co) the 
subtle competition between kinetic and Coulomb ener- 
gies together with the Pauli principle leads to the oc- 
currence of a ferromagnetic phase. Since the transition 
happens when both contributions to the total energy are 
of the same order, theoretical methods to study such 
systems have to be nonperturbative/ Within the dy- 
namical mean-field theory (DMFT);2*i2*Ai which is gen- 
uine nonperturbative, detail conditions for occurrence 
of ferromagnetism in a one-band Hubbard model were 
elucidatediiSiiii The same nonperturbative scheme to- 
gether with a density functional theory in the local den- 
sity approximation was used to describe ferromagnetic 
phases of Fe and NiiAi The presence of a binary alloy dis- 
order introduces a new nonperturbative aspect into the 
problem. Namely, when the difference between ion ener- 
gies A is larger then the band-width W, the conducting 
band is split and two alloy subbands are formed. 1 - 

Recent investigation of the one-band Hubbard model 
with a binary alloy disorder demonstrated an intrigu- 
ing interplay between effects due to interaction and 
randomness^ It was shown in Ref. [lfj that for elec- 
tron densities n < x, where x was a fixed concentra- 
tion of an alloy ion, and for large local Coulomb inter- 
actions U the Curie temperature T c is enhanced when A 
is increased. Additionally, at n — x the Mott-Hubbard 
metal-insulator transition (MIT) occurs for A = A C (U). 
The Mott-Hubbard transition in a non-integer filled sys- 
tem with binary alloy disorder was investigated further 
in Ref. where the notions of an alloy Mott insulator 



(A > U) and an alloy charge-transfer insulator (A < U) 
were introduced. 

In the present paper we study the one-band Hub- 
bard model and determine its physical properties, such 
as the Curie temperature, as a function of x. The Mott- 
Hubbard MIT can be reached by varying x while n, A 
and U are constants. Such MIT will be called an alloy 
concentration controlled Mott transition. This notion 
extends the previous classification of band-width con- 
trolled and filling controlled Mott transitions^ As will be 
presented below, under certain circumstances the Curie 
temperature has a maximum at a finite alloy concentra- 
tion. Additionally, new quantum phase transitions be- 
tween disorder ferromagnetic and disorder paramagnetic 
ground states in the vicinity of the Mott insulators are 
identified^ These quantum phase transitions are driven 
by varying the alloy concentration while the electron den- 
sity is constant. All those effects, as explained below, are 
caused by a subtle interplay between alloy band splitting 
and correlations. 

The present work is motivated by real experimental 
situations where the tuning parameter is rather x than 
A. Namely, while the latter is fixed for given atoms the 
former can be varied by making different alloy composi- 
tions or changing chemical stoichiometry. Thereby, the 
controlling of Curie temperature T c (x) is experimentally 
accessible. The examples are bec Fe-Co or fee Ni-Cu 
alloysi 19 : 20 In the first case, due to alloying, the system 
is driven from weak to strong ferromagnetism, whereas in 
the latter case the system is changed from a ferromagnet 
to a paramagnet. Interestingly, the behavior of the sat- 
urated magnetization and the Curie temperature in Fe- 
Co are non-monotonous functions of x reaching maxima 
at 30% concentration of Co>iS Another interesting ex- 
amples are pyrite alloys Ti_ 2; T' 2 ;S2, where T (T') stands 
for transition metals: Fe, Co, Ni, Cu, ZniSi Starting from 
FeS2 and alloying it with Co, the empty e g band (FeS2 
is a narrow-band semiconductor) is progressively filled in 
with electrons and the system becomes a disordered itin- 
erant ferromagnet. Surprisingly, the maximum T c (x) in 
Fei_^Co^S2 occurs not at x = 1 but at x ss 0.76. In dif- 
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ferent compound UCu2Si2- a: Ge a: the Curie temperature 
has a maximum at x w 1.6, i.e. again when the system is 
disordered^ Interesting aspect of this alloy is that since 
Si and Ge are isovalent, the system is isoelectronic and 
the disorder has only structural character not influencing 
directly the magnetic arrangement. 

In Section II we introduce a one-band Anderson- 
Hubbard model with quenched binary alloy disorder and 
solve it within DMFT framework. In Section III we 
present our results on ferromagnetic properties of this 
model when the density of electrons is independent of al- 
loy concentration, i.e. the system is isoelectronic when 
x is varied. Next, in Section IV we discuss our results 
when the system is non-isoelectronic under the change 
of x. Section V presents conclusions and a final discus- 
sion. 



II. MODEL AND DYNAMICAL MEAN-FIELD 
THEORY 

A. Anderson— Hubbard Hamiltonian 

Itinerant electron ferromagnetism in binary alloy is de- 
scribed hereby the Anderson-Hubbard Hamiltonian with 
uncorrelated on-site disorder 



H 



ein i(T + %^»1 
i 



(1) 



where is the hopping matrix element, U is the local 
Coulomb interaction, cf a is the fermionic creation op- 
erator for an electron with spin a in Wannier state i, 
and fiicr is the particle number operator. The quenched 
disorder is represented by the atomic energies e,, which 
are random variables. We consider binary alloy disorder 
where the atomic energy is distributed according to the 
binomial probability density 

P(e i )=x5(e i +jj+{l-x)sL-jj. (2) 

Here A is the energy difference between the two atomic 
energies, providing a measure of the disorder strength, 



while x and 1 — x are the concentrations of the two 
alloy atoms. For A 3> W, where W is the band-width, 
it is known that binary alloy disorder causes a band 
splittingii 5 ^ 2 ^ The number of states in each alloy subband 
is equal to 2xN a and 2(1 — x)N a , respectively, where N a 
is the number of lattice sites and a factor two counts the 
spin degeneracy. 



B. Dynamical mean— field theory 

The Anderson-Hubbard Hamiltonian is not solv- 
able at any finite space dimension. It is however nu- 
merically solvable in an infinite dimension after a proper 
rescaling of the hopping parameters^ i.e. a set of self- 
consistent DMFT equations is derived Aiiii The local 
nature of the theory implies that short-range order in 
position space is missing. However, dynamical correla- 
tions due to the local interaction and disorde r 25 ! 26 are 
fully taken into account. 

In the DMFT scheme the local Green function G an is 
given by the bare density of states N°(e) and the local 
self-energy Y* an as 



G an = / de 



N°(e) 



lUJ n 



(3) 



Here the subscript n refers to the Matsubara fre- 
quency u> n — (2n + 1)tt/P for the temperature T, with 
j3 = 1/ksT, and /i is the chemical potential. Within 
DMFT the local Green function G an is determined self- 
consistently by 

r I $ D[c <Jl c a ]c an c an e^K^} \ 

aU= \ fD[c„,c*]e^^} / dis ' W 

together with the k-integrated Dyson equation 

San = Gan + ^c™- (5) 

The single-site action Ai for a site with the ionic energy 
€i = ±A/2 for i =A and B, respectively, has the form 



dTC* a (T)Ca(T)c*_ a (T)C-a(T), 



(6) 



where we used a mixed time/frequency convention for 
Grassmann variables c CT , c* . In the presence of binary al- 
loy disorder, the single impurity problem has to be solved 
twice in each self-consistency loop. Averages over the 



randomness are obtained by 26 

<•■•>*. = / deP(e)(---), (7) 



which is equivalent to treating a problem of disorder 
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within the coherent potential approximation scheme^^ 
Due to the local nature of the theory and the arithmetic 
averaging of the physical one-particle quantities, Ander- 
son localization is not capturedi^LSS 

An asymmetric density of states is known to stabilize 
ferromagnetism in the one-band Hubbard model for mod- 
erate values of Jj£J2J£ Therefore, we use the density of 
states of the fcc-lattice in infinite dimension, 29 



N°(e) 



exp[_ibJs 



tt(1 + V2e) 



(8) 



This density of states has a square root singularity at 
e = — l/v2 and vanishes exponentially for e — > oo. In 
the following, the second moment of the density of states 
is used as the energy scale and is normalized to unity. 
Since the sa me density of states has been used earlier in 
Refs. |l2Hf| we are able to compare the specific numerical 
results. 

The one-particle Green function in Eq. J2J is deter- 
mined by solving the DMFT equations iterativelyiSii^ 
using Quantum Monte-Carlo simulations with the Trot- 
ter slice At = l/4i2Si Since we are mostly interested in 
qualitative behavior of T c vs. x at different U, n, and 
A we do not perform the extrapolation of the results to 
At — > 0. Curie temperatures are obtained by the diver- 
gence of the homogeneous magnetic susceptibility explic- 
itly implying that the ferromagnetic phase transitions are 
of the second order 



III. RESULTS FOR ISOELECTRONIC ALLOYS 
A. Curie temperature 

The Curie temperature as a function of alloy concen- 
tration exhibits very rich and interesting behavior as is 
documented in Figs. Q and [21 It is usually expected that 
T c (x) is suppressed in disordered systems. This indeed 
is found for most cases when x is varied between zero 
and one. However, at some concentrations and certain 
values of [/, A and n, the Curie temperature is enhanced 
above the corresponding value for the non-disordered 
case (x — or 1). This is shown in both panels of Fig. 2] 
for 0.4 < x < 0.9 and in the upper panel of Fig. [21 for 
< x < 0.2. The relative increase of T c can be as large 
as 25%, as is found for x w 0.1 at n = 0.7, U = 2 and 
A = 4 (upper panel of Fig. EJ) ■ 

This unusual enhancement of T c is caused by three 
distinct features of interacting electrons in the presence 
of binary alloy disorder^ 

i) The Curie temperature in the non-disordered case 
T c p = T C (A = 0), depends non-monotonically on band 
filling n. Namely, T^(n) has a maximum at some filling 
ii = n*(U), which increases as U is increased^ see also 
our schematic plots in Fig. [3J 

ii) As was described above, in the alloy disordered sys- 
tem the band is split when A ^> W . As a consequence, 
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FIG. 1: Curie temperature as a function of alloy concentration 
x at U = 2 (upper panel) and 6 (lower panel) for n — 0.3 and 
disorder A = 1 (dashed lines) and 4 (solid lines). 
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FIG. 2: Curie temperature as a function of alloy concentration 
x at U = 2 (upper panel) and 6 (lower panel) for n — 0.7 and 
disorder A = 1 (dashed lines) and 4 (solid lines). 



for n < 2x and T -C A electrons occupy only the lower 
alloy subband and for n > 2x both the lower and upper 
alloy subbands are filled. In the former case the upper 
subband is empty while in the later case the lower sub- 
band is completely full. Effectively, one can therefore 
describe this system by a Hubbard model mapped onto 
the either lower or the upper alloy subband, respectively. 
The second subband plays a passive role. Hence, the sit- 
uation corresponds to a single band with the effective 
filling n ff — n/x for n < 2x and n e s — (n — 2x)/(l — x) 
for n > 2x. It is then possible to determine T c from the 
phase diagram of the Hubbard model without disorder. 12 

iii) The disorder leads to a reduction of T^(n e g) by 




FIG. 3: Schematic plots explaining the filling dependence of 
T c for interacting electrons with strong binary alloy disorder. 
Curves represent T p , the Curie temperature for the pure sys- 
tem, as a function of filling n at two different interactions 
U\ <€. U21— Upper panel: For n < x, T c of the disordered 
system can be obtained by transforming the open (for Ui) 
and the filled (for U2) point from n to n e fi = n/x, and then 
multiplying T p (n/a;) by x as indicated by arrows. One finds 
T c {n) < T p (n) for Ui, but T c (n) > T p (n) for Ui. Lower panel: 
For n > x, T c of the disordered system can be obtained by 
transforming T^{n) from n to n e fl = (n — 2x)/(l — a;), and 
then multiplying T p [(n — 2a;) /(l — x)x] by 1 — x as indicated by 
arrows. One finds T c (n) > T p (n) for Ui, but T c {n) < T c p (n) 
for U 2 . 

a factor a = x if the Fermi level is in the lower alloy 
subband or a = 1 — x if it is in the upper alloy subband, 
i. e. we find 

T c (n) « «T c p (n cff ), (9) 

when A 3> FT" (C.f. Appendix). Hence, as illustrated in 
Fig- El T c can be determined by T^(n e s). Surprisingly, 
then, it follows that, for suitable U and n Curie temper- 
ature of a disordered system can be higher than that of 
the corresponding non-disordered system [cf. Fig. E] . 

The explanation of the T c (x) enhancement, given 
above, is supported by a detailed analysis how T c changes 
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FIG. 4: Changes of the Curie temperature with disorder A 
at x — 0.1 and n = 0.7 for different interactions U. 

when A increases at fixed x. The numerical results are 
shown in Fig. 0] at x = 0.1 and n — 0.7; ex amp les at 
x = 0.5 already have been presented in Ref. 16]. For 
n > 2x (results in Fig. ^corresponds to this regime) the 
Curie temperature initially decreases upon increasing A 
from zero. However, when A > U the trend is inverted 
and T c increases, finally saturating. At A ~ U the al- 
loy band splitting becomes effective, changing the be- 
havior of T c versus A. As shown in Fig.0] only for small 
U the Curie temperature is elevated above the value at 
A = 0. This is strongly related to the non-monotonic 
dependence of Tg(n), and in particular to the fact that 
its maximum changes with U. Namely, as is illustrated 
in the lower panel of Fig. El only at small interactions 
T% (n c ff) > TP(n), which is a necessary condition for the 
enhancement of T c by disorder. In the n < 2x case, on the 
other hand, the necessary condition T^(n e s) > TJ"(n) for 
an enhancement of T c implies that the interaction must 
be strongii& 

B. Magnetization and Curie constant 

The two different cases n < 2x and n > 2x (when 
A ^> W) should correspond to two different behav- 
iors of a saturated magnetization, i.e. the magnetiza- 
tion for T — > 0. Namely, in the first case only the 
lower alloy subband is occupied and the magnetization 
M = (n-j-) — (hi) — n, where (n^) are the average num- 
bers of electrons with spin up or down. However, in the 
second case M = n — 2x since the lower alloy subband 
is split off and fully occupied by the electrons with two 
spin species, thereby being magnetically neutral. 

In order to confirm this physical picture we calculate 
the magnetization M(T) as a function of temperature. 
The resultant magnetizations together with the inverse 
static susceptibilities X" 1 ^) are presented in Fig. [S] for 
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FIG. 5: Magnetization and inverse static susceptibility at: 
x — 0.6, n = 0.3 and U = 6 (circles and dashed lines), and 
x = 0.1, n = 0.7 and U = 2 (squares and solid lines). In both 
cases A = 4. 

two cases: i) x = 0.6, n — 0.3, and U = 6 (filled 
and open circles), corresponding to n < 2x case, and 
ii) x — 0.1, n — 0.7, and U = 2 (filled and open 
squares) corresponding to n > 2x instance. At both 
parameter sets T c is enhanced by the alloy band split- 
ting (A — 4). The numerically calculated magnetization 
(filled circles and squares) follow very closely the theo- 
retical Brillouin curves (dashed and solid lines) ~ The 
magnetization data, shown in Fig. are consistent with 
our conjecture that the saturated magnetization should 
be: M = n = 0.3 in case (i), whereas M — n — 2x = 0.5 
in case (ii). 

Two interesting observations are made: Firstly, in the 
case (ii) (with n > 2x) , the presence of disorder increases 
T c while the saturated magnetization is suppressed below 
its value at x — 0. In other words, the disordered sys- 
tem becomes a weaker ferromagnet but with higher T c . 
Secondly, although the correlated electrons in the disor- 
dered system are itinerant, the magnetization M(T) is 
well reproduced by the Brillouin curve, albeit formally 
this curve is derived for localized moments^ The last 
observation calls for an analytical proof within DMFT 
analogous to that already given for the linear behavior of 
the inverse susceptibility^! 

Within DMFT the uniform spin susceptibility for a 
systems with paramagnet-ferromagnet phase transition 
obeys the Curie-Weiss law x(T) = C/(T - T C )£L where 
C is a Curie constant. The same Curie- Weiss law is 
also derived with a mean-field theory for localized mag- 
netic moments and in this case the Curie constant is di- 
rectly proportional to the saturated magnetization. To 
check whether similar relation holds for itinerant elec- 
tron system within DMFT we compute the Curie con- 
stants obtaining: C\ = 0.240 ± 0.002 in the case (i), and 
C*2 = 0.385 ± 0.001 in the case (ii). It turns out that the 



ratio C1/C2 = 0.623 ± 0.005 is very close to the ratio of 
saturated magnetizations 0.3/0.5 = 0.6. This result pro- 
vides a new interpretation of the Curie constant within 
DMFT, i.e. C is proportional to the saturated magneti- 
zation, similarly as in localized magnetic moment theory. 
In addition this finding corroborates our initial conjec- 
ture regarding saturated magnetization in two physically 
different cases n < 2x and n > 2x at large A. 

C. Metal— insulator transition 

Upon increasing U and A, the Mott-Hubbard MIT 
occurs at the electronic filling n = xi&H Such MIT 
can also be encountered by varying x. This MIT will be 
called an alloy concentration controlled Mott transition. 
Approaching a correlated insulator the itinerant ferro- 
magnetism is also suppressed due to localization of the 
electrons. Using linearized DMFT 33 it was estimated^ 
that when A — > 00 the critical interaction U c « 6y/x. 
Such estimation leads to U c « 3.3 for n — 0.3 (Fig. QJ, 
and to U c s» 5 for n = 0.7 (Fig. [2J- This means that at 
A = 4 and U — 6 (lower panels in Figs. ^ and [2} the 
Mott-Hubbard MIT is possible at x = n. 

D. Quantum phase transitions 

At large U and A the system encounters several quan- 
tum (T = 0) phase transitions when x is varied. One is 
the Mott-Hubbard MIT discussed above. The other is 
the transition from a ferromagnetic metal to a paramag- 
netic metal. In this quantum phase transition the Curie 
temperature vanishes at x = x c which is different from 
x = n, where Mott-Hubbard MIT occurs. 

The presence of the quantum phase transitions should 
affect the properties of the disordered Fermi liquid at fi- 
nite temperatures, provided that the system is of finite 
dimensionality. This scenario would be worth to be in- 
vestigated including spatial correlations, which however 
goes beyond the present DMFT framework. 

IV. NON-ISOELECTRONIC ALLOY 

In many alloys the change in the concentration x al- 
ters the electron density n. In this Section we investigate 
T c (x) when x and n are varied simultaneously. The re- 
sults presented in Fig. are obtained under the assump- 
tion that x = 2n. When x = the band is empty (the 
system is a band insulator) . Upon increasing x from zero 
to one, the band filling increases from zero to one-half 
(quarter filled band). This model realization would cor- 
respond to the physical situation in Fei_ a; Co a ;S2 alloy if 
the number of states in e g band is normalized to one. Of 
course, this analogy should not be stressed too far since 
in our model an important exchange interaction (Hund's 
coupling) is absent. 
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FIG. 6: Behavior of the Curie temperature for U — 2 (upper 
panel) and 6 (lower panel) when the change of electron filling 
n is associated with the change in the alloy concentration x. 
At large U the Curie temperature is enhanced in the disorder 
system. 

As shown in Fig. the presence of disorder at weak 
interaction always suppresses the Curie temperature with 
respect to the non-disorder case. However, at large U 
and A > 1 the Curie temperature is larger than that in 
the pure case at the same filling. 

This difference again can be understood on the basis 
provided by the scheme depicted in Fig. |3J In the 
present case n < 2x for all x, and only the lower alloy 
subband plays a role in the effective description at large 
A. In this limit this subband is effectively filled with 
ri e ff = n/x = 1/2 electrons. Using the same arguments 
as in the upper panel in Fig. [3] we see that T c (x) is 
enhanced only for large U . 



V. CONCLUSIONS 

In the present paper we studied the one-band 
Anderson-Hubbard model with binary alloy disorder 
showing that the Curie temperature in such alloyed corre- 
lated electron system can reach higher values than those 
in the non-disorder system. We also identified and dis- 
cussed the metal-insulator transition at non-integer fill- 
ings. Additionally possible quantum phase transitions 
at zero temperature were described pointing out on a 



possible interesting quantum critical behavior of a disor- 
dered Fermi liquid. This work completes previous studies 
of itinerant ferromagnetism in the pure and disordered 
Hubbard model within DMFTi 12 i 13 i 16 

Regarding the physical systems, it is now of great 
importance to extend the one-band Hubbard model 
with alloy disorder to a multi-band case and formulate 
the adequate version of DMFT for solving this problem. 
It is very interesting which aspect of the ferromagnetism 
in one-band alloy are generic and will be present in 
multi-band case. 
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APPENDIX: HARTREE-FOCK THEORY 

Within the Hartree-Fock approximation one can find 
analytically that in the strong disordered limit A ^> W 
the self-energy has the form 

UM , X A a;(l-a;)A 2 , K . 

En. = CT— + (1-2*)-+ aUM V* A » (^.1) 

4 Z z 4 

where M is the magnetization density. This self-energy 
leads to the splitting of the density of states with x and 
1 — x of the initial states in each alloy subband. Since 
the integrated function in the equation for T^ F is peaked 
at w = /i, only one of the subband gives contribution 
to evaluate T C HF . As a result, T C HF rj qT c hf p, where 
a = x or 1 — x, as introduced in Section II. Although 
in DMFT one cannot find analytically the correspond- 
ing self-energy, the splitting of the density of states also 
appears and only one of the subbands contributes in the 
equation for T c . In analogy to the Hartree-Fock approx- 
imation, we assume that T c is reduced by a with respect 
to T^, which we find to be valid even at strong interac- 
tion. 
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